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1 Introduction 



High energy e~^e~ annihilation experiments have proven to be a particularly clean testing 
ground for quantum chromodynamics (QCD), the theory of strong interactions. One of the 
key predictions of QCD has been the occurence of jets of hadrons in high energy reactions Jl], 

Quantitatively, the observed hadronic events can be classified according to the number of 
jets. The most commonly used definitions of jets in e~^e~ annihilation are based on iterative 
jet clustering algorithms. One begins with a set of final-state particles, and considers for 
each iteration the pair {i,j} with the smallest value of a dimensionless measure yij. If this 
Uij is smaller than a preset cutoff i/cut, the pair is clustered into a single pseudoparticle 
according to a recombination rule. The procedure is repeated until yij > ycut for all pairs of 
(pseudo-)particles, at which point the remaining objects are called jets. In the perturbative 
calculation, the quantity y^j is constructed from the parton momenta. 

There exist a number of jet observables that are well defined in QCD, and which can 
be calculated perturbatively as an expansion in the strong coupling a^. It is also possible 
to define a variety of quantities (so-called event shape variables) that depend on the final 
state structure and do not need the above concept of a jet. One of the well-known event 
shape variables is thrust The only requirement for such quantities is that they are 
infrared safe, i.e. they neither distinguish between a final state with one particle being soft 
and a final state where this particle is absent, nor between a final state with two coUinear 
partons and a final state where these two collinear partons are clustered into a single hard 
pseudoparticle. Infrared safe quantities that are dominated by three jet configurations, and 
the three jet differential cross section (defined by the jet algorithm described above) are 
of interest, in particular because these quantities are in leading order (LO) proportional to 
as- The next-to-leading order (NLO) QCD corrections to the production of three partons 
were computed more than fifteen years ago |§, |, ||, 0] for massless quarks. Subsequently, 
these results have been implemented in numerical programs [§-|TB[ and widely used for 
tests of QCD with jet physics. Also, a complete calculation for oriented three jet events 
with massless quarks was performed in p3[ . 

To date huge samples of jet events produced at the Z resonance have been collected 
both at LEP and SLC. From these data large numbers of jet events involving b quarks can 
be isolated with high purity using vertex detectors. While the mass of the b quark can 
be neglected to good approximation in the computation of the total bb production rate at 
the Z peak, the n-jet cross sections depend, apart from the cm. energy i/s, on an extra 
scale, namely the jet resolution parameter ycut- For small ycut effects due to quark masses 
can be enhanced. It is therefore desirable for precision tests of QCD to use NLO matrix 
elements that include the full quark mass dependence. One can then also try to extract 
the mass of the b quark from three jet rates involving b quarks at the Z peak, as suggested 
in [|16|, elaborated in [|1^], [Q, and experimentally pursued by the DELPHI collaboration 
ITOj . Further applications include precision tests of the asymptotic freedom property of 
QCD by means of three jet rates and event shape variables measured at various center-of- 



mass energies, also far below the Z resonance For theoretical predictions concerning 
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the production of top quark pairs at a future e+e~ collider, the inclusion of the full mass 
dependence is of course mandatory. 

The three, four, and five jet rates involving massive quarks have been computed to 

In reference 



leading order in already some time ago [21, 22 



23 , the NLO corrections to 



the production of a heavy quark pair plus a hard photon are given. Recently, results for the 
decay rate of the Z boson and the e^e~ annihilation cross section at the Z peak into three 
jets including mass effects at NLO have been reported Jl^, ||2^, HTSf. Further, momentum 
correlations in Z ^ bbX at NLO have been studied [p5| ]. 

We have computed the complete differential distributions for e+e^ annihilation into three 
and four partons via a virtual photon or Z boson at order a^, including the full quark mass 
dependence. This allows for order predictions of oriented three jet events, and of any 
quantity that gets contributions only from three and four jet configurations. In this article, 
we would like to discuss this calculation in detailQ. We start in section ^ with an overview 
of the calculation. Section ^ contains a general kinematical analysis and the LO results 
for e~^e~ — > 3 jets including mass effects. The virtual corrections are discussed in section 
|[ They contain singularities that cancel against real soft and collinear contributions from 
four parton final states. To isolate the latter contributions, a modification of the phase 
space slicing method incorporating massive partons is developed in section ^. In section 
^ we show numerical results. In particular, we compare for some observables the massive 
with the massless result. We summarize our results in section |^. An Appendix contains the 
analytical result for the virtual corrections to the three parton production rate. 



2 Outline of the calculation 

The calculation of an arbitrary quantity dominated by three jet configurations and involving 
a massive quark-antiquark pair to order proceeds as follows: 

We first have to compute the fully differential cross section for the partonic reaction 

e+e- ^-f*,Z* ^QQg (2.1) 

at leading and next-to-leading order in a^. Here Q denotes a massive quark and g a gluon. 

At leading order, the parton momenta can be identified with the jet momenta. If the 
identity of the particles making up the jets is known, we have for the differential three jet 
cross section 

dcri(e+e" ^ 3 jets) = 0(?/qq - ycut)'d{yQg - ycut)'d{yQg - ycut)d(Ti{e'^e~ QQg), (2.2) 

where the lower index 1 means the result in order a J, 6 is the Heaviside step function, ycut 
is the jet resolution parameter and yij G {yQQ,yQg,yQg} is determined by the experimental 
jet definition. For example 

yij = — ^(1 - cos%) for the JADE algorithm [29], 



short account of our work is given in 
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i/ij = ^(1 — cos^jj) for the Durham algorithm ||30||, (2-3) 



where 6ij is the angle between jet i and j in the cm. system. At higher order and in the 
experimental analysis, these definitions are supplemented by a recombination rule to form 
a pseudoparticle k from a clustered pair {i,j}- For the above algorithms, one simply adds 
the four- momenta, pk = Pi +Pj- (For a detailed discussion of these and other jet algorithms, 
see also If the particle content of one or more jets is not known, the LO calculation 



of a three-jet quantity from dai{e~^e~ QQg) may involve an averaging over the different 
ways of assigning jet momenta to parton momenta. 

The one-loop integrals, which appear in the NLO amplitude for ( ^.11 ), contain ultraviolet 
(UV) and infrared (soft and collinear) (IR) singularities. In the framework of dimensional 
regularization in = 4 — 2e space-time dimensions, these singularities appear as poles in e. 
The UV singularities are removed by renormalization. Details will be given in section ^. 

After renormalization, the virtual corrections to the differential cross section for (p.l|) still 
contain IR singularities. These have to be cancelled by the singularities that are obtained 
upon phase space integration of the squared tree amplitudes for the production of four 
partons. The relevant processes are 

e+e" 7*,Z* -> QQgg, 
e+e" 7*, Z* QQqq, 

e+e- ^^*,Z* ^QQQQ, (2.4) 



where q denotes a light (massless) quark. The first two reactions of (|2.4| ) give singular 
contributions in the three jet region when a parton becomes soft and/or two partons become 
collinear. The cancellation of the singularities has of course to be performed analytically. 
To this end, the soft and collinear poles of the relevant four parton matrix elements have 
to be explicitly isolated. We do this by modifying the so-called phase space slicing method 
lOi to account for masses. (For alternative methods, see, e.g., 0, |12|, |l^ and references 



therein.) 

The idea of the phase space slicing method is to introduce an unphysical parameter Smin, 
which is much smaller than all relevant physical scales of the problem, e.g. Smin ^ si/cut 
for jet cross sections defined by the jet resolution i/cut- The parameter Smin splits the four 
parton phase space into a region where all partons are "resolved" and a region where at 
least one parton is soft and/or two partons are collinear. In the massless case, a convenient 
definition of the resolved region is given by the requirement Sij > Smin for all invariants 
Sij = {pi + VjY- We will modify this definition to account for masses, see section ||, but 
still will use the terminology "resolved" and "unresolved" partons. In the regions with 
unresolved partons, soft and collinear approximations of the matrix elements, which hold 
exactly in the limit Smin 0, are used. The necessary integrations over the soft and collinear 
regions of phase space can then be carried out analytically in d dimensions. One can thus 
isolate all the poles in e and perform the cancellation of the IR singularities between real 
and virtual contributions, after which one takes the limit e — > (cf. eq. ( |5.25| )). The result 
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is a finite differential cross section (i(T^(e+e" — > QQg) at order a1 for tliree resolved partons, 
wliicli depends on Smin- Its contribution to a given tliree-jet quantity is specified by tlie 
experimental jet definition, wliicfi may involve an energy ordering, a tagging prescription, 
etc. 

The contribution to a three jet quantity or to an event shape variable of the resolved 
part of the four parton cross sections (i(T^(e^e^ 4 partons) is finite by itself and may 
be evaluated in c? = 4 dimensions, which greatly simplifies the algebra. Its contribution 
to three jet configurations also depends on Smin and is obtained by recombining the four 
resolved partons into three jets according to some recombination scheme and by projecting 
the four parton phase space onto the experimentally specified three jet phase space. 

Since the parameter Smin is introduced in the theoretical calculation for technical reasons 
only and unrelated to any physical quantity, the sum of all contributions to any three jet 
observable must not depend on Smin- In the soft and collinear approximations one neglects 
terms which vanish as Smin ~^ 0. This limit can be carried out numerically. Since the indi- 
vidual contributions depend logarithmically on Smini it is a nontrivial test of the calculation 
to demonstrate that the sum of resolved and unresolved contributions becomes independent 
of Smin for small values of this parameter. Moreover, in order to avoid large numerical 
cancellations, one should determine the largest value of Smin which has this property. 

Schematically, we have for the differential three jet cross section involving massive quarks 
at order a^: 

rf(T2(e+e- 3 jets) = 0|rfa^(e+e- QQg) 

da^ ie^e^ QQgg) + (icr^(e^e~ QQqq) 

+ daf (e+e- QQQQ) |, (2.5) 

where contains the experimental jet definition, and the integration symbol contains 
the procedure of recombining resolved partons into jets and of projecting the four parton 
phase space onto the three jet phase space. Eq. ( p^.5D is schematic in the sense that the 
exact relation between a given three jet quantity and the parton cross sections depends on 
the experimental knowledge about the jets, as already mentioned above. 

It should be emphasized that the notion of resolved/unresolved partons is unrelated to 
the physical jet resolution criterium or to any other relevant physical scale. In particular, 
we will use the freedom in defining the soft and collinear regions of phase space to simplify 
the necessary analytic integrations in d dimensions. Explicit expressions for the resolved 
cross sections entering ( |2.5| ) will be derived in section ^. 

The fully differential cross section for the production of two jets with massive quarks is 
known to order as already for some time For order a1 calculations of quantities that 



also get contributions from two jet configurations, the two loop amplitude for e"'"e~ — QQ 
is needed. Results for the total cross section for e^e~ hadrons to order with mass 
corrections are reviewed in ||2^. (For recent results see also [^.) Since 

a{e^e~ — > hadrons) = a{e^e~ — > 2 jets) + cr(e^e^ 3 jets) 
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+ a(e+e" jets) + 0(<), (2.6) 

one can now also calculate the two jet cross section to this order, but not yet arbitrary 
differential distributions. 



3 Kinematics and leading order results 

In this section we review some of the basic kinematics relevant for the process 

e+(p+)e-(p_) Z*,7* - Q{k,)Q{h)g{h) (3.1) 

and give the leading order results for the fully differential cross section for ( p. 11 ). The 
amplitude T for the reaction ( |3.1| ) (see also Fig. [^) can be written, to order af^ and to 
leading order in the electroweak couplings in the general form: 

+ v{p+)^,u{p^){-QqV^)]. (3.2) 



In ( p.2| ), s = (p+ + P-Y, Qq denotes the electric charge of the quark in units of e = v47ra, 
and (7^, gl are the vector- and the axial-vector couplings of a fermion of type /, i.e. 

gl = Ti -2Qfsm^^w, and gl = Ti . 

In particular, g*^ = — | + 2 sin^ ^w, g'a = ~\ for electron or (/^ = — ^ + | sin^ ^dw, ga = ~\ 
for a bottom quark, with dw denoting the weak mixing angle. The function x(s) is given 

by 

1 s 
^^^^ ^ 4 sin^ §w cos2 s-ml + imzTz ' ^^'^^ 

where rrtz and Tz stand for the mass and the width of the Z-boson. We have defined in 
( p.2| ) the amplitudes V^, A^, and , which contain the information on the decay of the 
vector boson into the three final state partons. The contribution is generated at order 
a^J"^ by a quark triangle loop The sum runs over the different quark flavors. Only 

quark doublets with a mass-splitting give a non-vanishing contribution. At higher order in 
as there are additional terms in eq. ( p.2|) of the form V^" ^ which are absent at order a^J"^ due 
to a non-abelian generalization of Furry's theorem. Such additional terms are also present 
in the scattering amplitude for the four parton final state qqQQ with light quarks at the 
initial vertex at lowest order in a.^- 

The squared amplitude |Tp can be written to order in the compact form: 



iri 



s2 



LPO^u^pc ^ L^^'^'^fff/l . (3.4) 
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We restrict ourselves here to the relevant case of longitudinal polarization of electrons and/or 
positrons and neglect the lepton masses. The lepton tensors L^*^^" and L^^^'^ take the usual 
form: 

LPC,u ^ = pIp':_ +pIp>^ - g'^'^p^p., 

LPV,.^ = ^IV[75j^+7^;^_7.] = -i£'^>^p^, (3.5) 

where £0123 = +1- The tensors H^^^^^^ (which we call "parton tensors") may be written 
as follows: 



ttPC(PV) _ A^V ttVV j_ AA ttAA , VA+ ttVA+ , AA^ ttj 

^tiu — 9pC{PV)^ij.u ^ 9pC{PV)^tJ,v ^ 9pC{PV)^nv ^ 2^9pc{pv)^h 

i 

, V- TT^^X , VA- rrVA. , H-^^- (o(i\ 

+ / .9pn(pv)^iJ-^ ^9pc{PV)^nv + / .9pn{pv)J^iJ-^ ■ l^^-^J 

i i 

The coupling constants are given explicitly by 

9pc{pv) = Qq fpc{pv) + 2 Qq R'ex(s) fpc{pv) + fl'?^ lx('S)|^ fpc{PV)^ 

9pC{PV)—9a \X\S)\ JPC{PV)^ 

9pi(pv) = -9a Qq R'ex(s) fpc{PV) - 9v 9a \x{s)\'^fpc{PV), 

,AAi O „i |,,/„m2 nZZ 



9pC{PV) — 9a 9a\X\o)\ JPC{PV)^ 

9pciPV) = -9a Qq Re x{s) fpl\PV) - 9v 9^ \x{s)\'^fpS{PV), 
9pi{pv) = ^9aQQ Imx(s)/pc(py)' 



FA' 



9pc{PV) ^i9lQQ Imx(s)/pc(py), (3-7) 



where 



fpc — 1 ~ ^-^+^ 

f/^ = -(l-X.X^)g: + {X.-X^)g:, 

f^S = (1 - A-A+)(^f + gf) - 2(A_ - A+)^X, 

JPV — ^- — 

r/y = (1 - X^X^)g: - (A_ - X^)g:, 

/l^ = (A- - K){gf + gf) - 2 (1 - A„A+)^?: g^ (3.8) 

with A_ (A_|_) denoting the longitudinal polarization of the electron (positron) beam. 

The parton tensors , H^^, . . . arc related to the amplitudes V^, and in the 
following way (we do not write down explicitly the sum over the polarizations of the final 
states) : 

= 



6 



AA 



jjVAi 
rrVA± 



V,Al±A^V:, 



A 4" 



flU 

In massless QCD chiral symmetry implies that if^^ 



4* A* 



vv 



H^^^ and i/;/- 



(3.9) 

0. The tensors 



H^jf- and Hl^u' multiply factors that are formally of higher order in the electroweak 



.VA 



coupling. Their contributions to the differential cross section for (|3.1| ) are negligibly small. 

The parton tensors H^^ entering ( p. 41 ) contain the complete information about the orien- 
tation of the partons in the final state with respect to the beam axis and about their energy 
distribution. Although the parton tensors are good for calculational purposes^, we find it 
more convenient for the phenomenological application to switch to a slightly different nota- 
tion, which we adopt from ref. We write the e~^e~ —>■'-/*, Z* —>■ QQg fully differential 
cross section (in d = 4 dimensions), with no transverse beam polarization, as follows: 



3 Q!o 



dx dx d cos 6 ( 



^pt 



47r TT 
-|- F4 sin 26 cos 
-|- F7 sin 29 sin 



Fi(l + cos^ ^) + ^2(1 - 3 cos^ e) + F3 cos 



) + -F5 sin^ 6 cos 2(f) + Fq sin 9 cos < 
-|- Fg sin^ sin 20 + Fg sin 6 sin cf 



where 



a 



pt 



cr e e 
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47ra' 



(3.10) 



(3.11) 



In ( p.lO|) , X = 2kki/s and x = 2kk2/s are the scaled energies of the quark and antiquark 
in the Z rest system, 6 is the angle between the electron direction and the quark direction, 
and is the (signed) angle between the e~^e~Q plane and the QQg plane. The functions Fi 
are easily obtained by contracting the parton tensors H^^, with suitable tensors built from 
the final state momenta and ^f^j^. 

The functions F7, Fg and Fg are generated by absorptive parts in the scattering ampli- 
tude, i.e. they are zero in tree-approximation. In addition they vanish for massless quarks 
at one loop order [Q, ||14||P|. For the nonvanishing LO functions F^ we have: 



F' 



(0) 
.(0) 



^VVp(0)yV AAp{0),AA 



1,2,4, 



5), 



VAj, 

9pv 



(J = 3,6). 



Using the abbreviations 



m 



s 



B 



1 



(1 -x)(l -x)' 



2 — X — X, 



(3.12) 



(3.13) 



^For the computation of iJ^jy in the massless case to order a1, see 



^ A nonzero function gives rise to a correlation between the beam axis and the normal to the three jet 
event plane. This "event handedness" was first discussed by and studied in detail at the Z resonance 
both experimentally ISa and theoretically p6|. 
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where m is the quark mass, we find: 



(0),VV 



B 



+ z 



3(x + xy + 8(x + x) + 2xx{l -Xg)-Q B - 2z^XgB L 



Fl 



(0),AA _ TpiO),VV 



n{o),vy 



^'^^ + + x)^ - 10(1 - Xg) + QzxlBy 

^-^|l - Xg- z{-3x^ + x^ - 2xx + 2x-2x + 2)B 



x"^ — 



2z^[2x^ - 13x^ - + 2xx(4x + 3x - 11) + 20x + 12x - 8]5^ + Sz^^x^B'^ k 



3^2 d2 



,(0),AA 



1 - Xg z{l - x) [x^ - - + xx{2x + x - 8) + 8x + 6x - 2] - 22;^x^ 



X^ — 4:Z 



;i-x)2(x2-42) 



2B 



VX^ — 4^ 



x^ + 2x^ + xx(-x + 2) - 2x - 22fx^ - 8x^ - 4x^ 



+ xx(6x + 5x - 16) + 12x + 8x - A]B + Sz^x^B \, 



Ao),vv 



BJ{1 - x)il - x){l - Xg) - zxl ( 



(x2 -4^)(x2 -4z) 



X [2(1 — Xg) — xx] 



2z 



-[Ax +xx(3x + XX — 6x + 6) — 4x — 8x + 4 

1 — X 

4z(x^ + x^ — 2x^ + xx(x + 3) — 4x) — IGz'^Xg] 



p(0),AA ^(o)yv , 4zi?^(l-x)(l-x)(l-xg)-2:xg 

Fl " = H <x -2(l-x„)+xx 

x^ — 4^; ' 

^io),vv ^ ^io),vv _ 2^^!^ ^^2^ 



^(0),VA+ _ jl - X){1 - X){1 - Xg) - ZX I J ^ 

^ V x2 - 4^ 



X 



1 — X 



(3.14) 



Our results agree with the resuhs of ||3^ and reduce in the massless case to those given in 
p2|Pi The increase in algebraic complexity due to the non-neglection of the quark mass is 
quite substantial already at LO. 



4 Virtual corrections 



The calculation of the virtual corrections for reaction (|3.1| ) is straightforward albeit tedious. 
Non-neglection of the quark mass leads to a considerable complication of the algebra. In 

* In reference [Q, the expression for F^^'''^^ — p^^')'^^ (\a,st equation in (Bl) of denoted there 
pio),a _ p{o),as^ contains a misprint: The denominator of this function is given there as (1 — a;)^(l — x)^, 



while the correct denominator is (1 — x){l — x), cf. eq. (3.14) above 
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this section we give details on this part of the calculation. 

We work in renormalized perturbation theory, i.e. start with a renormalized Lagrangian 
including counterterms to calculate the (truncated) renormalized Green function for reaction 
(|3.1| ) in d dimensions to order a^/^. The relevant one- loop diagrams are shown in figure 0, 
the counterterm diagrams are depicted in figure S. The scattering amplitude, 



e e 



QQg) = T 



Born, 



e e 



virtual / 



e e 



(4.1) 



follows from the renormalized truncated Green function by multiplication of UV finite wave 
function renormalization factors in accordance with the Lehmann-Symanzik-Zimmermann 
reduction formalism. Except for the quark mass renormalization, we work in the modified 
minimal subtraction (MS) scheme throughout. For the quark mass renormalization we 
consider both the on-shell scheme and the MS scheme. The relation between both schemes 
is given to order by the well-known relation 



pole 



MS 



a^(/i) / 4 



2 In 



TT 



m^) 



(4.2) 



where [i is the renormalization scale. It is somewhat more transparent to use the pole 
mass from the start and, if desired, switch to the MS scheme only at the very end of the 
calculation using ([4.21). In an ah initio calculation with the MS mass one has to distinguish 
between the renormalized mass parameter of the propagators (which then by definition is 
the MS mass) and the mass of the external particles (which of course is always the pole 
mass) . 

Throughout the calculation presented in this work the background field gauge pSI, p9| is 



used. This gauge simplifies the triple-gluon vertex and in addition a simple relation between 
the renormalization constants of the coupling Zg and the wave function of the gluon field 
Za holds in the MS and MS schemes, namely 



-1/2 



(4.3) 



This implies also that the renormalization constants of the quark-gluon coupling Z^p and 
the quark wave function Z^ are equal. 



(4.4) 



For the calculation of the loop diagrams the Passarino-Veltman method |4^ is used to 
reduce the tensor integrals to scalar one-loop integrals. When doing the loop integration 
UV as well as IR singularities are present. Dimensional regularisation is employed to treat 
both. The formulas of [H0| are generalized in order to account for the appearance of IR 



poles. The reduction of tensor integrals and the trace algebra in d dimensions is carried out 



with two different programs: FoRM [41| and REDUCE |42| are used independently for this 
part of the calculation and yield the same results. 
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Because of the presence of the axial- vector current, a prescription to handle the 75 
matrices in d dimensions has to be chosen. In the case of traces with two 75's, we work with 



an anticommuting 75 in d dimensions, which is known to be consistent By using the 

relation 7I = 1 we eliminate the 75's from the traces. In the case of traces with only one 
75 the situation is more complicated. Here we use the 't Hooft-Veltman prescription |^ 



It is well known that this prescription violates certain Ward identities, i.e. in the limit of 
vanishing quark masses chiral invariance is broken. To restore the chiral Ward identities, a 
special 75 counterterm has to be taken into account in higher order calculations. Explicitly, 
in our case the replacement [4^ 



7 (y 
7^75 - ^5'TTTe;..,.,.37^^7'^V^ with = 1 - (4.5) 

O! TT 

restores the chiral Ward identities to order a^^^. 

When calculating the loop integrals we find it useful to replace the box integrals in 
0? = 4 — 2e dimensions by the box integrals in 6 — 2e dimensions plus a linear combination 



of the three-point integrals. This method |46| has the advantage that the IR singularities 



only appear in the three-point integrals, because the box integrals in 6 dimensions are IR 
finite. In this way the IR divergent part of the virtual corrections can be easily obtained. 
Further, the coefficient functions which multiply the loop integrals become algebraically less 
involved in this way. 

We compute both the real and the imaginary part of the loop integrals by Feynman 
parametrization techniques. The imaginary parts of the integrals induce several contribu- 
tions to the differential cross section ( |3.10|) : The function Fg is nonzero only through these 



imaginary parts, and the parity violating functions F-j and Fg get the dominant contribu- 
tions from the imaginary parts of the one-loop integrals^. Our results for the imaginary 
parts of the integrals agree with the results given in ||36[] . 

The explicit expressions for the virtual corrections to Fi — Fg are too lengthy to be fully 
reproduced hereQ We will restrict our discussion of the analytic results to the function Fi, 
which determines in particular the three parton rate, 

(jptFi. (4.6) 



dxdx 71 

The virtual corrections to Fi may be written as 



ipvirtual VV ipvirtuaiyV . AA Tpvirtual,AA . -^r;^ AAi Tpvivtual,AAi / A ^\ 

^1 =9pc^i +9pc^i +Z^9pc^i ■ 14. g 



^Additional contributions to i^7,8j which are formally of higher order in the electroweak coupling, are 
induced by the imaginary part of the Z propagator times the real parts of the one-loop integrals via 7 — Z 
interference. While the main virtual contributions to the functions ^3 and Fg are proportional to the real 
parts of the one-loop integrals, there are additional subdominant contributions from j — Z interference of 
the form (imaginary part of one-loop integrals) x (imaginary part of the Z propagator) . 

^The result for Fg is given in [p6[. 
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We can further decompose the first two terms in ([4.7| ) as follows (recall that we already 
removed the UV poles by renormalization; the remaining poles are due to IR singularities): 



virtuniyv (AA) 



2nY 



2 1 
+ - 



e2 e 



17 



+ 2 In 



+ \n{zB) - 7 



?>N, 



c 



1 1 



(2/3-(l + /?2)lnM 



(0),VV{AA) 



j^rcst, singular, Vy{Ayl) j^counter, finite, yV(ylA) j^ext., finite, yy{Ayl) 
+ -^1 +-^1 +-^1 



+ 0{e 



(4. 



where 



/9 



I- Xa 



1+/? 



(4.9) 



n™*^ is the number of massless flavors, the functions F\ 
B and Xg have been defined in eq. ( p.l3| ). Further, 



1-/5' 

{Q)yv(AA) given in (|3.14|) , and 



rest, singular, yy 



la, n 



rest, singular, AA 



e Att 



NcxlBil + 2z). 



(4.10) 



The finite contributions from the counter terms and the external wave function factors are 
given in the on-shell mass renormalization scheme for one massive and n™*^ massless flavors 
by: 



counter, finite, yy 



a. 



B 



X 



„2 r:)/_3 



ln(47r)-7 + ^-lnfp 



x^ - zB{x'^ + Ix^ + 5xx(-x + 2) - 20s + 7) 
4z^E2(-4x^ + 16x2 + 2x^x2 + xx(3x2 - 17x + 17) - 23x + 6) 
Sz^XgB'^il - xf + {x^x] 



— Xg — ZXgB 



F counter, finite, AA y-icounter, finite, yy j~, 
1 ~ -'^1 

TT 



ZS 



ln(47r) - 7 + - - In — 
3 V/^ 



3(^2 - 8x + XX + 3) - 2zB{x^ - 2Ax^ + 3xx(x - 6) + 66x - 28) 



+ 2Az^XgB'^{l- xf + (x 



+ 3xo(x + x) -2zxlB 



Fl 



ext.,finite,yy 



as 
12^ 



llA/'c - 2n?')(ln(47r) - 7) - 21n 



ZS 



F^ 



(o),yy 
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lext., finite, 



(lliVc 
a, 



J I o 



zs 



-{nNc-2nT) 



llNc- 2nT) (ln(47r) - 7) - 2 In — 



{0),AA 



xlBil + 2z) 
f ' 2 



(4.11) 



The rather involved results for the finite leading-color [Fi 



lc,yy(AA)^ 



and subleading-color 



^p^c,vv{AA)^ contributions from the interference of the loop diagrams, Figs. 2(a)-(k), with 
the Born graphs are listed in the appendix. Finally, 



AAi T-ivirtual,y4Ai 



2 rZZ^s X -\- X 



PC 



TT 



[l-zBxl]Y.gl^^F[ 



jj, jj, 



(4.12) 



where Zi = mf / s denotes the square of the scaled mass of the quark in the fermion triangle of 
Figs. 2(1), (m) and Y.i fl'a-^e/*(x, x, Zi) can be obtained by simple substitutions from formulas 
(2.17) and (2.18) of ||3^ and will therefore not be given explicitly here. This contribution 
to Fi is finite by itself and numerically very small [p2| . 



5 The phase space slicing method for massive quarks 

In this section we will describe how we modify the phase space slicing method in the presence 
of masses and derive explicit expressions for the resolved cross sections rfo"^ at NLO entering 
equation (p75|). 

In the presence of massive quarks, the structure of collinear and soft poles of the four 
parton matrix elements is completely different as compared to the massless case. In par- 
ticular, the nonzero quark mass serves as a regularizer for collinear singularities. Thus, the 
matrix elements contain fewer singular structures, but the presence of a quark mass leads 
to more complicated phase space integrals. 

The contribution of the process e^e~ — > QQQQ to the three jet differential cross section 
is free of singularities. Thus, for this subprocess, it is possible to define: 

da^{e+e~ ^ QQQQ) = da^ie+e' ^ QQQQ). (5.1) 

The process e^e'^ QQqq is singular in the three jet region only when the massless partons 
become collinear, whereas the contribution of e^e~ —>■ QQgg contains both soft and collinear 
divergencies. We will first discuss the contributions of soft gluons. 

The amplitude for T{e~^e~ Q{ki)Q{k2)g{ks)g{k4)) can be written in terms of color- 
ordered subamplitudes: 

T(e+e- ^ QQgg) = (r'^«T"^),,e,^i + (r'^^T»^)c,c,^2, (5.2) 
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where denote color matrices, 03,04 the color of the gluons, and ci,C2 the color of the 
quarks. For the squared matrix element (summed over colors and spins) we thus get 



E 1^1 



QQ99W 



colors 
qqg-spins 



E 



qqg — spins 



N, 



c 



(5.3) 



The term 1 5*1 + 5*21^ in (|5.3| ) contains the QED-like contributions. 

In the limit where one gluon becomes soft, each of the terms in the squared matrix 
element ( |5.3|) can be written as a factor multiplying the squared Born matrix element for 
e+e' QQg: 



Nl-IN, 



c 



Born , 



e e 



qqg — spins 



colors 
qqg-spins 



fc4-o glNc 



/K2,4,3) \r''°'\e^e- ^QQg)\\ (5.4) 



colors 
qqg-spins 



Here Qs denotes the strong coupling constant. The limiting behavior of the term |S'2p follows 
from ( |5.4| ) by the exchange (1 ^ 2). For the term subleading in the number of colors we 
have 



Nl-l 
4iVc 



E \Si + s, 



qqg — spina 



1 2 fc3-*0 



fc4^0 



91 



2Nc 

,2 



/.,(!, 3, 2) Y l^^°™(e+e- 



colors 
qqg-spins 



9 s 



2N, 



c 



colors 
qqg-spins 



We have defined in (|5.4| ), { f>.5\ } the eikonal factors 

Atik 4m2 



t 



2 ' 



4ti2 4m2 4m2 



/2 ' 

t9„ 



//(^,s,^) 
fsiil,s,2) 

i^lsi^2s 

where tij = 2kikj and m denotes the quark mass. 
For notational convenience, we further define 

Qisfc = 0('^js ~l~ tsk '2^min)y 
^isk — 0(2s,njn tis tsk)- 

Bearing in mind the soft limits ( |5.4| ), let us consider the identity 

i^ir = (6134+6134) (6243+6243) 

= (61346243 + 6134 + 6243 ~ 61346243) I'S'ip. 



QQ9)\', 

QQg)\%5.5) 



(5.6) 



(5.7) 



(5.^ 
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The first term in the second line of eq. ( [5 .81) corresponds to the case where the energies 
of both gluons are bounded from below, i.e., neither of the gluons can become soft. (They 
can still become colhnear, which will be discussed later.) The second and the third term 
give rise to soft gluon singularities in the three jet region. Finally, the last term describes 
the situation where both gluons are soft. This term therefore does not contribute to the 
differential three jet cross section, but to the C(a^) two parton cross section. (Also the 
second and the third term contain such contributions and the negative sign of the fourth 
term compensates the double counting.) 

We will now derive the complete contribution from soft gluons to the three jet cross 
section at NLO. We start by discussing the contribution from [S'lp. With our choice of 
the Heaviside functions defining the soft region, it is simple to analytically integrate the 
eikonal factors ( |5.6| ) in d dimensions over the soft gluon momentum. Let us consider the 
case where ^3 ^ (second term in ( p.8|) ). The case —>■ can be treated in complete 
analogy. Since the two gluons are identical particles, the overall soft factor multiplying the 
Born cross section for e~^e~ — > QQg is determined by considering the contributing from one 
soft gluon only and leaving out the identical particle factor (1/2!). In the cm. system of 
the heavy quark Q{ki) and the hard gluon g{k4), the eikonal factor reads 



//(I, 3, 4) 



1 



EiEi 




2Vti4 + m'^ 



m 



cos ^^)(1 — /?i4 cos 6') -Ei(l — /3i4 cos I 



(5.9) 



where Ei^^ are the heavy quark and soft gluon energies in that system, 9 is the angle between 
the heavy quark and the soft gluon, and /5i4 = ti4/ (ti4 + 2m?). In the same system we have 
(with (i = 4 - 2e) 



134 



e( 



^7r2 r(l - e 



- E^-^' {sin 9)-^' dEsd cos 9. (5.10) 



The integration over the soft gluon momentum k^ can now be carried out without difficulty. 

The soft contribution from 15*2^ is obtained in the same manner by switching the roles 
of the heavy quark Q{ki) and the heavy antiquark Q(^2)- For the term subleading in color, 
we write 



\Sl + S2\ 



(ei32 + 6132) (ei42 + 6142) |5i + S2\^ 
©1320142 + ©132 + ©142 — ©132©142 ) I'S'i + 5*21^. 



(5.11) 



Since the amplitude Si + 5*2 is a QED-like contribution with massive quarks, it does not 
induce colhnear but only soft singularities |^ . For ^3.4 — we choose the quark- antiquark 
cm. system to evaluate the eikonal factors. The complete soft factor S multiplying the 
squared Born matrix element =°i"s |T^°'''^(e^e^ — > QQg)\^ which we obtain by adding 



colors 
qqg — spins 



all contributions reads, if we relabel the remaining hard gluon momentum with ^3: 



S = ^Nc ^ 

4:71 r(l - e) I Smin 



ti3 + 



In 1 + 



t 



13 



14 



+ 2 ln(2) - 1 

- iln= 

2 



TT^ r 

— + 2 \n\2) - 2 ln(2) + 2 ln(2) + ^ + 1 

6 ^^13 



ln(l + ^^ 



2Li2 



13 



ti3 + m ^ 
2/3 ' 



13 



t 



23 J 



+ 2 ln(2)(l + p^) In(cu) - In^cu) - 2(1 + /52)Li2 



1 + ln(o;)J - 4/3 ln(2) + 2 In(cu) 

■ 2(3 \ 



1 + 



+ 0(e), (5.12) 



where (3 and a; have been defined in terms of the scaled quark and antiquark cm. energies 
X, X in ( |4.9| ). (We have ti2 = s(l — Xg — 2z), tis = s(l — x), t23 = s(l — x).) 

The leading color contribution l^ip + 15*21^ to the matrix element (|5.3|) also contains 
coUinear singularities. We isolate them by writing 



01340243|'S'l|^ — [©1340243 " 0(i^^l3 ~ 2Smm)0(i'^24 — 2Smm)0(Smm 
+ 0(^13 ~ 2Smj„)0(t24 ~ 2Smm)0(Smm ~ ^34)|'S'lp. 



';34j 



l^ll 



(5.13) 



The first term in ( [5.13| ) is now free of singularities over the whole four parton phase space 
allowed by the Heaviside functions, while the second term contains the contribution of two 
coUinear gluons which are both not soft. By construction we thus avoid an overlap of the 
soft and the coUinear part of phase space. The resolved part of the differential cross section 
for e~^e~ QQgg entering (|2.5|) may thus be defined as: 



1 1 



da^ie^e ^ QQgg) = — - 



n 



X 



2s 2! [fj^ i2Tif2Ei_ 
Nl-lNc 



{2Ti)%k-Y.k, 



i=l 



I — -T72"01320142|5'i + 5*2 

5 — spins ^ C 



+ [01340243 ~ 0(^13 ~ 2Smm)0(^24 ~ 2Smj„)0(s^ 
+ [02340143 - 0(^23 " 2S™n)0(tl4 - 2S™n)0(Sr 



^34) 
^34) 



l^l 
152 



(5.14) 



The remaining calculation is completely analogous to the massless case [IDj; we include 
it here for completeness. 



In the limit fc-? II ^4 we define 



with kh = + k^. In this limit, 

INc ^ , ^ ,0 fc3||fc4 gsNc 



Born / 



e e 



(5.15) 



QQg)\\ (5.16) 



qqg- spins 



colors 
qqg — spins 
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where 



/ 



99^9 



2 l + + 



(5.17) 



is proportional to the Altarelh-Parisi sphtting function. The colhnear behavior of |S'2p is 
identical to ( ^.16 ), leading to an overall factorization of the squared matrix element ( ^.3|) 
in the colhnear limit. Further, with tis = ^tih, ^24 = (1 — C,)t2h in this limit, 



0(^13 ~ 2Sm,j„) 0(^24 ~ 2Smm)0('Smm ~" ^ 



34 J 



(27r)'^-i2E3 (27r)^-i2E4 

0(^^1/1 — 2Smiri)0((l — 0^2/1 — 2Smm)0(Smm 

1 {AnY r. d'^-^kh 



t34) 



167r2r('l -e^ 



[t34e(l-e)]"'^^« 



{27iY-^2Eh 



(5.18) 



After integration over ^ and ^34, summing the contributions from l^ip and |5'2p, and 



relabelling kh ^ k^, we get for the colhnear factor C^^ multiplying ^oiois |T 
QQ5')P statistical factor 1/2! is included): 



Born I 



e ' e 



qqg— spins 



rv 

C99 = ^Nc 



27T 

3 



rri 



In(ei) + ln(6) + ^ 

D 



67 1 



1^ 



In2(e2) +0(6) 



(5.19) 



where ^1 = 2Smin/tl3, 6 = '2Smin/t23- 

As mentioned before, the contribution of the process e+e~ — > Q{ki)Q{k2)q{k3)q{ki) to 
the three jet cross section is singular only when the two massless quarks become colhnear. 
To isolate the singular term, we use 



|r(e+e- ^ QQgg)P = [0(^34 - s^in) + 0(s™„ - ^34)] |r(e+e- ^ QQqq)\^. 
For the resolved part of the cross section we may therefore write 

rfo-f (e+e" QQqq) = 0(t34 - s™„)(icr2(e+e~ QQqq)- 
In the phase space region defined by 0(smm — ^34), we use the colhnear limit 



(5.20) 



(5.21) 



e e 



QQqq)\ 



„2„ms 
-^|2 '=3l|fc4 {Js"-/ jqq^g 



Borii/ 



e e 



gg^?)p, (5.22) 



qqg -spins 



qqg— spins 



where 



/ 



2 ^2^(i_^)2_^ 

^34 1 - e 



(5.23) 
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and n™^ is the number of massless flavors. After integration over the coUinear phase space 
we get for the coUinear factor C"^^ multiplying the squared Born matrix element for e"*"e^ 
QQg: 



271 f 



'47r/i^ 



r(i 



(5.24) 



We have thus derived the following differential cross sections for three resolved partons 
entering (|2.5|): 



da. 



e e 



QQg) = TT 

2s 



n 



\ {2Tif2Ei_ 

X lim lr^°™ 



e-+0 I 



e e 



1=1 

. QQ(7)r"'-*"^'(e+e- QQg)* + h.c. 



+ 



S + C3^ + c^^ 



IT- 



Born I 



e e 



QQ9)\' 



(5.25) 



In particular, we have for the function Fi deflned in ( |3.1(]| ), ( |4.(j| ), which determines 
double differential cross section for three resolved partons: 



(1) 



with given in (CT ), ( CT ), and 



(1) 



lim 



virtual i Tjisoft+coUinear I 
+ ^1 f) 



where _pvirtuai j^^^g been listed in 



( PID and 



Fl 



soft+coUinear 



One can easily verify that F^"* is flnite (the poles in F™^^^"^ and Fl 



qq 



(0) 



■ yyxlB ^^ x]B{l + 2z) 
9 PC o 3 PC o 



(5.26) 



(5.27) 



(5.28) 



;isoft+collinear 



exactly cancel). 



The dependence of on Smin cancels in its contribution to an observable quantity (like 
the three jet cross section) against the s^m dependence from the contributions of four 
resolved partons in the limit Smin 0. The latter are obtained numerically starting from 
expressions ( p.lD , (|5.14| ), and ( ^.21| ). Examples will be given in the next section. We have 
derived analytic expressions also for all the other functions -F/^'', z = 2, . . . , 9 which enter the 
fully differential cross section (|3.10| ) for three resolved partons at NLO; but we will discuss 



them elsewhere [48 



6 Numerical results 

In this section we show results for some observables involving massive quarks. We carry out 



the necessary numerical integrations of our matrix elements with the help of Vegas 
All quantities are calculated by expanding in to NLO accuracy. 
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The three jet cross section for b quarks as a function of the jet resolution parameter i/cut 
in the JADE and Durham scheme as well as an observable sensitive to the mass of the b 
quark at the Z pole have already been presented in [|^. Here we start our discussion by 
demonstrating the independence of physical quantities on the parameter Smin as Smm 0. 
We choose as an example the three jet fraction for b quarks, 

tAycut) = — I — • (6.1) 

'-'tot 

In ( |6.1| ), the numerator a\ is defined as the three jet cross section for events in which at least 
two jets containing a 6 or 6 quark remain after the clustering procedure. This requirement 
ensures that the cross section stays finite also in the limit m;, — > 0. The contribution of 
the process e^e~ — >■ Z, 7* — >■ qqg* qqbb to the three jet cross section with one tagged b 
quark develops large logarithms In(m^) - which find no counterpart in the virtual corrections 
against which they can cancel - when the hb pair is clustered into a single jet. In principle, 
there are two distinct possibilities to handle this problem: One may either impose suitable 
experimental requirements/cuts to get rid of events with two light quark jets and one jet 
containing a bb pair (the definition for a\ chosen by us is an example for this), or one can 
improve the fixed order calculation by absorbing the large logarithm into a fragmentation 
function for a gluon into a heavy quark. A detailed discussion of this issue will be presented 
elsewhere [Q. Note that a^^^ has to be calculated only to order for the NLO prediction of 



/|; hence the prescription to handle the g* —>■ bb contribution does not affect the evaluation 



of the denominator of (|6.1| ) at this order. 

Figs. ^ and |^ show the three jet fraction /| in the JADE and Durham scheme at NLO 
as a function of i/min = Smin/{sycut) at a fixed value of i/cut = 0.03 and y/s = niz = 91.187 
GeV. The error bars are due to the numerical integration. For the renormalization scale 
we take in these plots /i = y/s. As to the mass parameter, we use m^^(/i) defined in the 
MS scheme at the scale /i. The asymptotic freedom property of QCD predicts that this 
mass parameter decreases when being evaluated at a higher scale. (A number of low energy 
determinations of the b quark mass have been made; see for instance Q and 

references therein.) With m^^(/i = mt,) = 4.36 GeV [|52| and as{mz) = 0.118 |^ as an 



input and employing the standard renormalization group evolution of the coupling and the 
quark masses, we use the value m^^(/i = mz) = 3 GeV. One clearly sees that /| reaches a 
plateau for small values of ymin- The error in the numerical integration becomes bigger as 
ymin ^ 0. In order to keep this error as small as possible without introducing a systematic 
error from using the soft and collinear approximations, we take in the following i/min = 10~^ 
for the JADE algorithm and ymin = 0.5 x 10~^ for the Durham algorithm. At these values, 
the dominant Smm-dependent individual contributions from three and four resolved partons 
are about a factor of 2.5 (JADE) and 4 (Durham) larger than the sum. 

In Figs. I and |^ we plot /| as a function of ycut at LO and NLO, again at y/s = mz- 
The QCD corrections to the LO result are quite sizable as known also in the massless case. 
The renormalization scale dependence (where /i is varied between and 2m^), which 

is also shown in Figs. ^ and ^ is modest in the whole ycut range exhibited for the Durham 
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and above Ucut ~ 0.01 for the JADE algorithm. Below this value perturbation theory does 
not yield reliable results in the JADE scheme. In Figs. ^ and ^ we take a closer look on the 
scale dependence of /|, now using the on-shell mass renormalization scheme. We vary the 
scale 11 between m^/lO and 2mz for a fixed value Ucut = 0.2 x 10~^/^° ~ 0.1 and on-shell 
masses = 3 GeV and m^"^'^ = 5 GeV. In Fig. ^ we see that the scale dependence 

of the LO result (which is solely due to the scale dependence of as at this order) in the 
JADE algorithm amounts to about 100% in the fi interval shown. The inclusion of the 
corrections reduces the scale dependence significantly; the NLO result for /| at /i = 2m^ 
is about 30% smaller than the NLO result at /i = m^/lG. In the case of the Durham 
algorithm, the difference between /| at /i = 2mz and at /i = m^/lG is reduced from about 
100% at LO to about 10% at NLO. 

The effect of the b quark mass may be illustrated by looking at the double ratio 

CiVcu.) = J^^y (6-2) 

where the denominator is the three jet fraction when summing over all active quark flavors, 
which is given to a very good approximation by the massless NLO result @]- Similar 
double ratios have been studied in and |T6|, Ijl^, ||T8[. In Fig. |10|we plot C as a function 



of the cm. energy at Ucut = 0.08 for the JADE algorithm. The running of as is taken into 
account in the curves, where we again use as an input as{fi = mz) = 0.118. For the dashed 
(LO) and full (NLO) curve we use the running mass with m^^(/i = mz) = 3 GeV. For 
all energies, the renormalization scale is set to fi = a/s. For comparison we also show the 
LO result for a fixed value of the b quark mass = 4.7 GeV (dash-dotted curve), which 
is the corresponding value of the pole mass. One clearly sees that the effect of the b quark 
mass gets larger for smaller cm. energies. 

Another interesting quantity to study mass effects is the differential two jet rate J56| 
defined as 

n f ^ /2(l/) - f2{y- Ay) 

D2{y) = , (6.3) 

where f2{y) is the two jet fraction aX y = Ucut for a given jet algorithm. The advantage of D2 
over the three jet fraction /s lies in the fact that the statistical errors in bins of D2{y) are 
independent from each other since each event enters the distribution only once. To order 
a^, D2{y) can be calculated from the three- and four jet fractions using the identity 

l = /2 + /3 + /4 + 0(a3). (6.4) 

We define 

D'2iy) 
DT\y)' 

where we - as in the case of the quantity C - use the massless NLO result to evaluate the 



ny) = (6-5) 



denominator. We plot our results for in Figs. Oand O, again for i/i = /i = mz- The 
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full circles show the NLO results for m^^(/i = mz) = 3 GeV. For the O(a^) contribution of 
fiiUcut) to V we use a fit to the numerical results. For comparison, the squares (triangles) 
are the LO results for rrib = 3 GeV (m^ = 5 GeV). The horizontal bars show the size of the 
bins in ycut- The effects of the b quark mass are of the order of 5% or larger at small values 

of Vcut- 

Finally we discuss an event shape variable which does not depend on the jet finding 
algorithm, namely thrust , defined as the sum of the lengths of the longitudinal momenta 
of the final state particles relative to the axis n chosen to maximize this sum, 

T = max " L (6.6) 

We may write the thrust distribution for b quarks as 



dT 271 ' \2n 



In the massless calculation at /i = ^/s, the coefficients Ci and C2 are independent of the 



cm. energy. This changes in the massive case as shown in Figs. ^ and [Tj. For these and 
the following plots of the thrust distribution we exclude the singular two-jet region near 
T = 1. We use a slicing parameter s^m = 1 GeV^. We now only require the tagging of 



one b quark and omit the contributions from the process e~^e~ — > qqg* qqbb. In Fig. [15 
we plot the thrust distribution ( |6.7|) at ^/s = /i = 30 GeV with ^^(/i = ^/s) = 0.142 and 
T^b^^il^ = y/s) = 3.33 GeV. Shown separately (and not included in the NLO histogram) 
is the contribution from e~^e~ — > qqg* —>■ qqbb calculated "naively", i.e. directly from the 
matrix element without imposing cuts. The LO and NLO thrust distributions for b quarks 



at a/s = mz depicted in Fig. |T6| are almost identical to the corresponding distributions at 



a/s = 30 GeV. The reason is that the change in the coefficients ci^2 due to mass effects as 



exhibited in Figs. |13| and |Tj is almost exactly compensated by the larger value of as at 



fi = 30 GeV. The contribution from g* bb splitting is bigger at the higher energy scale. 



7 Summary 

In view of the large number of jet events collected both at LEP and SLC it is desirable for 
precision tests of QCD to use NLO partonic matrix elements that include the full quark mass 
dependence. We have presented in this article the necessary ingredients to calculate any 
jet quantity that gets contributions from three- and four jet final states involving massive 
quarks at order a^. In particular, we have derived an explicit analytic expression for the 
virtual corrections to e~^e~ QQg- Our approach involves a modification of the phase space 
slicing method incorporating massive quarks. The independence of physical quantities on 
the slicing parameter s^m is a crucial test of the overall consistency and has been carefully 
checked. As to numerical results we have evaluated the three jet fraction for b quarks and 
compared it to the inclusive three jet fraction evaluated for massless quarks at different cm. 
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energies. We have further studied the dependence of our results on the renormahzation scale. 
If we take this dependence as an estimate of the neglected higher order corrections, we may 
conclude that they are of moderate size, in particular in the case of jet rates computed 
using the Durham algorithm. A particularly interesting quantity to study mass effects is 
the differential two jet rate. At the Z pole and for small ycut^ the effects of the h quark mass 
are of the order of 5% or larger. As a final example we studied the thrust distribution for h 
quark samples at ^/s = rriz and y/s = 30 GeV. Future work will include the calculation of 
several other three jet quantities, in particular also parity violating observables. 
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A Appendix 

In ( f4.8| ), we defined the finite contributions p^'^'^^^'^'^^ and ^j^e function 

pvntua.i,vv(AA} ^ which are generated by the interference of the loop diagrams Figs. 2(a)-(k) 
with the Born graphs. In this appendix, we give explicit expressions for these contributions 
in terms of finite parts of one-loop integrals times coefficient functions. For the two- and 



three-point loop integrals we use the notation of |^0|, i.e. all integrals are labelled with 
respect to the following box integrals: 

° Z7r2 J {P + ie){{l + hY + ie){{l + k^f - m? + ie){{l - /ta)^ - m? + ie) ' 

^sc,i ^ ±_ f {2nfi)-''d'-'H 

° Z7r2 J {P + ie){{l + hf - m2 + ie){{l + ki^Y - nP + ie){{l - k^Y - rrP + le) ' 

^o^'' = ^o^''l.,^.„ (A.1) 

with kl = kl = m?-, fc| = 0, and ^13 = fci + ^3. In the following, Co(l,2,3) denotes the 
integral obtained from Dq by omitting the fourth denominator, etc. The nonvanishing one- 
point function is defined by A{zs) = {2'7ifi)~'^'' Jd'^l[P — + As stated in section ||, 
we eliminate the above scalar box integrals in ci = 4 — 2e dimensions in favor of the box 
integrals in ci = 6 — 2e dimensions. 
We write 

+ cr^^^^^^Re Co(l, 2, 3) + cr.I^^'^^ReColl, 2, 4) 
+ c'-r^^^^Re Co(l, 3, 4) + 4-r''^^''^Re Co(2, 3, 4) 
+ 6lf ^(^^)Re So(l, 3) + }}-^''^^^^ReM2, 4) 
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+ 6lf ^(^^)Re 5o(l, 4) + fer'^^^-'^Re 5o(3, 4) 



(A.2) 



sc,VV{AA) _ JVV{AA)j^^ j^sc,l,d=fi , ^,^^(^^) 



d=6 



ReDn 



+ d 



'Re a 



+ c^^^^^Re era, 2, 3) + ~cl,:'->Re 0^(1, 3, 4) 
+ cZ^^^^ReC-'\l, 2, 4) + d-r^^^^Re Co(l, 3, 4) 



~VV{AA). 



sc,2,d=6 


■ysc,2 / 



sc,VV{AA) 
234 



Re Co(2, 3, 4) + cta^^^^Re 0^(2, 3, 4) 



■ysc,l / 



+ c; 

+ 6r/''^^^^Re 5o(l, 3) + blf'^^^^ReBo{2, 4) 
+ 6r/''^^^^Re5o(l,4) + blf''^^^^ReBo{3,4) 



^24 



ReF;''(2, 4) + a^^'^^(^^)Re A + 



(A.3) 



The symbol / means that the finite part of the UV or IR divergent integral / is to be 
taken. 

The real parts of the integrals appearing in ( [A.2| ) and ( |A.3| ) read (recall that z = 
w? / s, Xg = 2 — X — x): 



(1 — x)a 



+ 



sXg \ {a — x%){a — xi) 
(1 — x)x\ vr^ 
(a — x%){x% — xi) 6 

' x\ — p\ ( 1 — p — X 



TT^ fa -I 

-— + Li2 

o \ a 



X , 



- Li, I 



1] +Li ( ^+ ~ ^ 



+ Lia ( rd_^ ) _ Li, 



+ 



1 — X 



Lig 



+ - In 



+ / 

1 — Xj 



1 — x^_ — 



2 yx:^ — 

(1 — x)xi 
(a — xi)(xi — x+) 
^xi — 



-In'. 

2 V 1 - p - x^ y 



x: 



+ Li2 



xl — p 



x: 



x: 



ReCo(l,2,3) 



+ Li2 

+ -ln2 

2 yxi — p/ 

+ (x ^ x), 

1 1 f 1 1 



Li, 



1 — p — X* 
1 — xi 



Lio 



-X 



1 — xi — p^ 



2 I 1 - p - x^ 



2s 1 - X e2 



7 + In 



'sz(C-l)^ 



47rp2 

-l('l57r2 + 24Li2 ('^') +121n2(C) 



(A.4) 
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6 



7 + In I 



\^ 47r//2 

ReCo(l,2,4) = ReCo(l,2,3) 



1 2 ■ 



ReCo(l,3,4) = — 



Xj 



Li2 



— p 



+ Li2 



\x% — p ^ 

x\ 

x% + p-l, 
'x^ + p — l' 



2 

* +Li2 



— 1 
^x%+ p-l^ 
'xl - 1\ 



Xj 



ReCo(2,3,4) = ReCo(l,3,4) 
1 



ReSo(l,3) = 



7 — In [ — ^— r I + ^— — In , ^ , , ^ 



^ V^ln(C) + 2, 



Re5o(2,4) = Re5o(l,3) 

1 

e 
1 



ReSo(l,4) = - - 7-ln 
Re So (3, 4) 



l^4V. 
/ zs 



Re A{zs) = zs ( 7 — Inf 



+ 2, 

+ 2 + (1 - 2p) In 



ReCr'(l,2,4) 



1 1 



1 (1,,/ .) 



si — 1 — 27]] e \1 — r] 



-2U2 



-In 



1 — T}^ 



lln 

2 



_ V 



V 



ReCr''(l,2,3) = - 



1 — 
1 1 



ReC^'''(l,3,4) = ReC^'''(l,2,3) 
1 fl 



ReC^^'^(2,3,4) 



SX g I 2 



In" 



P 



TT 



1 — 7^ 



2[{l-Xg){l-x){l-x)-zxl] 
X |(l-x)(l-xJ(l-X5-4z) 



X -x)ReDo' 



sc,l,(i=4— 2e 



+ ReCr(l,2,4)] 



(A.5) 
(A.6) 



(A.7) 
(A.8) 
(A.9) 
(A.IO) 
(A.11) 

(A.12) 

(A. 13) 



(A. 14) 

(A.15) 
(A.16) 
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where 



+ [{l-x){l-xg)-2zxg\ 

X [(1 - x)Re Co''^(l, 2, 3) - XgRe Co''^(2, 3, 4)] 

+ {l-x)[-x{l-Xg) +2z(x + x)]ReCo(l,3,4) 



ReL); 



sc,2,d= 





SC,1 / 



Ren 



sc,l,d=6 




ReS^^'^(2,4) = i - 7 - In + 2 + (1 - 2r^) In 



1 — ?7 



a 



c 



1 — x 

1 — X + z 



1 



p = -(i-v^r^), 



-(X ± Vx"^ — 4:Z), 

1 



1 - -{x ± Va;^ - 4z), 



'1 



1~ Xr. 



(A. 17) 
(A.18) 

(A.19) 

(A.20) 

(A.21) 
(A.22) 
(A.23) 
(A.24) 

(A.25) 



and (cf. |5 



Ren 



sc,l,d=4-2e 




1 1 



1 



1 



1 — X I — Xg 1 — 2r] 



In 



1 — 1] 



— TT 



- 2 Li2 \l 
+ Li2 (1- 



Tl P 

l-r]l- 



- 2 Lio 1 



V 1 - P 

l-T] p ^ 



(A.26) 



The coefficient functions read (where B = 1/(1 — x)/{l — x), and /3 and u have been 
defined in (O)): 



"d=6 



TT 

„2 



+ 4;z5 -2x^ + 5x + xx{x - 2) - 2 - z(2x'' + 3(1 - 2x) + xx) 



+ (x ^ x) 
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-123 



Y^ja;^ + 2zB[-?>x^ + + 2xx{x - 2) - 3 - zx^] + {x ^ 



'-124 ~ ''123 



ic,yy 

Cl34 



zs (-4:z{l + 2z) 2x'^ + 2xx - 12x + 9 - 4:X - x'^ 



1 — X { 1 — X 1 — X 

z[13x^ - 34a; - 12x + 8xx + 25 - 8^(2a; + x - 3)] 
(1 - x)^ 

+ — f— — x^x + 8x — 6xx + 9x — 7 

x^ -Az 

+ z{36x - - 60 + 29x - lOxx) + Sz^{x + x - 6)] 
3 

- ^ ^[-(1 - + z(8x^ - 20x - 16x + 20a;x - 5x^x + 12) 

{x^ - 4z)^ 

+ Az^{x'^ - lOx - 5x + 2xx + 12) + 16/] I , 

ic,yy _ ic,yy 

^234 — Ci34 

^ic.yy _ 3^ ^(1-a;) -I + 2z{b - 2x) + Az'^{x - A) + lQz^ 
~ l-x ~ 2(1 -x)2 ^ 2(1 -S)(l -4^) 

-l + a; + 2;2(l -x) +;22(^a;-2) + 2/ 1 



2(1 - X + z){l - z){l - x) 2(1 - a;)(l - z){l - Az){x^ - Az) 

X [-(1 - a;)(l -Xg + 3xx) + z{Z7 - 26x + lAxx - llx^x - 2x^ - 18x) 
+ /(124a; + 4a;^x - 112 + a;^ + lOOx - 48a;x) 

+ 2/(-120 + 2x^x + 5x + 10a;x) + 8/(2a;x + 32 - 12a; - llx) + 64/] 
TiZ 

- 1 ^[(1 - 2;)(2a; - 3a;x + 5x - 2) 

+ z{-6x^ - 12xx + x^x + 2Sx + 16x - 24) + 4/(-6 + 2a; + x)], 



u^VV _ f\c,VV 
^24 — ^13 

B 



±3 ( 

^icyv ^ hx^ -3x^ - x^ - AOx + xx(2x^ + 2x^ + 5xx - A3x + 36) + 8 

16 I 

+ 8zB[-x^ + 7x^ + 8x'^ - A8x + xx{x^ + x^x - Ax^ - 3xx - a; + 20) + 20] 
+ <oAz^B\-3x^ + 12a;^ - 18a; + xx{x^x + a;^ - 8a; - a;x + 11) + 5] 

+ 128/S%(1 -a;)^! 

xB 

+ ^—^ —\x^ - 8x^ + 2Ax^ - 30a;2 + 22a; - 6 

8(a;2 - Az) 

+ x[x^ + 4a;^ - 27a;^ + 48^^ - 44a; + 12) + x^{Ax^ - \7x^ + 22a; - 6)] 

3x\x^ ^ XX -2{\- Xg)f 
+ 16(1 -x)(x2- 4^)2 +(^^^)' 
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i/^^^ = —\2x^ + xx{2x^ + X - 6) - Sx^ - 2x^ - 12x + 8 
8 I 

+ %zB{xx{x^ ^xx- llx + 16) - + Vlx^ - 28x + 10 + 2zx\ 
^ - [2x^ - 2x^ - + 28x^ - 36^2 + 28x - 8 



x"^ — Az 

+ 2x{x^ + 2x^ - 18x2 + 32^2 - 30x + 8) + x'^{<ox^ - 2hx^ + 32x - 8)] | 
?>x\x^ ^xx-2{\-Xg)f 



8(1 -4z) 



2 



+ (a; a;), 



^ic.yy _ — i5[-llx^ + 24x + 2xx(4x-7)-7]+4zS2[4x2-2xV-16x^ + 23x 

+ xx^-Zx" + ITx - 17) - 6 - 2^x„(l - xf\ - y^^'^^^-J^^-^^^^ 
^ ' 9\ ) \ {2-x){x^-^z) 

2x^ - 8^2 + 12x - 7 + x{2 - xf 



2 ^ 



+ (x x) > , 

{2-x){l-x + z) ^ 'y 

H Ix"^ - I6x^ + 24a; + 3xx{x - 4) - 8 

8 I 



7 — m ' 



s(l-x)2 

+ 4^B[5x2 - Ibx'^ + 20x + xx{7x - 12) - 5] 

+ 16z'^B'^[-4:X^ + 2x^xp- + 16x^ - 23x + xx{?>x^ - Ylx + 17) + 6 

+ 2.x,(l - xf\ - ^fe!±^^|i^^| + _ (A.27) 



-123 



^AA 
<-124 



lc,AA 
Cl34 



2zSs[l - a:g - zBx\ 



TT 



a; 



10a; + a;a; + 5 



+ <azB{2x^ - 6a; + + 3) + (x ^ x)|, 



-123 

^AA 
-123 



2^5 

1 -a; 



X 



10a; + a;a; + 5 + 3zBXg + (a; 



a; 



^134 + 

1 — a; 



42;2 z(x'^ + x'^ - 16(x + x) + 2xx + 16) 



X 



3(1 - a;) 



-(1 - xY{l -x) + 2z{l - a;)(3a; + X - 9) + 2z'^{2x - x - 7) 



+ 



12z^{2x + x- 3) 



3(^2 - Az) 
2z 



[(l-a;)(a; + 3x-2a;x- 1) 



3(l-a;)2 (^2-4^)2' 
+ z{-hx^ + x^x + 20a; + 12x - 10a;x - 16) + 4z^(2a; + x - 5)] |, 



\c,AA \c,AA 

^234 — Ci34 
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ficAA _ ricyv / 1 + ^ ^ l-2z + Sz"" z(l - x) (1 - z)(l -x)-z^ 



1-x (l-x)(l-4z) {1-xY {l-x){l-x^z){\-z) 

1 r 9 

+ 2^(15x2 - 31a; + hxx - Ax + 13) + 2^2(-12x^ - 12a;x + llx + 28) 

+ Az^{Axx + 12x - 6x - 21) + 16^^] 
6 

+ 7 TT^T^ ^[(1 - + zi-Qx'^ + Sx^x + 16x - lAxx + 12x - 10) 

(1 — x)[x/ — Azy 



+ 2z^{-x^ + 12a; - 2xx + Qx- 16) - 8;z^]|, 



6; 



lc,AA 




24 




lc,AA 


ricyv 


14 


— 0i4 




— 14a;x - 



?>x'^[x'^ + xx-2{l-Xg)Y 
A I (1 -x)(x2 -4z)2 



+ %&z^B^x„il - xf - -^^—[^x^ - 19a;^ + 34a;^ - 30a;^ + 14a; - 2 

x^ —Az 

+ 2x{2x'^ - 3a; + 2)(2a;2 - 4a; + 1) - x^{l - x){2,x^ - 6a; + 2)] + (a; ^ x) >, 

l\c,AA i\cyV , 2 f 3.r[,r^ + ,7-,r — 2(1 — .Ty)]- ^d/o 2 on , o - , inN 

O34 =034 + o i rrr^ 25(3a; - 20a; + 3a;a; + 10) 



2 1 (1 -x)(a;2 -4z)2 

- SzB^xl + — [8a;^ - 29a;3 + 48a;^ - 36a; + 12 

^ x^ -Az 

+ 2x{6x^ - 19a;2 + 20a; - 8) + x'^{Ax'^ - 9a; + 4)] + (a; ^ x)|, 

1 r :r2 - .r - ,r.r + 2,r - 1 



{1 - X + z){2 - x){l - x) 



xBlx"^ + XX - 2{1 - Xg)]'^ _1 
+ 2(2-x)(a;2-4z) +(^^^)/' 



2(2 -x)(a;2 -4z) 



Attz/j, 



2 



+ 



[a;^ + a;a; — 2(1 — Xg)] 
x^ - Az 



2 



s(l-a;)2^ 

- lOa;^ + 40a; - 10a;x - 7 + 2zB{5x^ - 52a;^ + 15a;^x + 130a; - 46a;x - 52) 

- 2Az^B^Xg{l - xf + {x ^ x)y (A.28) 



4-6 = 7T—^ — ^-^^gJ \x2^2x''-2x-2x + 2xx + 1 

"•"^ 7r(l - x)(l - a;^ - Az) 
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0.7 

r 9 / \ 

[—3a; —Ax + 7a; + — 5 + xx[x + a; — 6) 



{l-x){l-xg) 





+ 4z(- 


~ vv 


— "d=6 


-yv 
•-123 


ZS 



yy _ ZS j 2(2 -x + z) 2z{l - x)[x{x - A) + Az{-3x + 2) + 16z^ 
1 — X (1 — a;)2(x — 4;i;)x 



(2a; - 1 - 4z)[(a; - 82)(2,x2 - 2a; + 1) + 8z'^{8x - 3) - 96^^! 



(1 - x)(l - X(, - Az){x - Az) 



2 



2 

+ 7 77 ^^[x^ix'^ + x + l) + 2zx^{2x^ - 4a; - 5) - 8^^a;(8a;^ - 3a; - 

{\ - x){x - AzYx^ 



134 



•^134 ~ "^123 



c[24 = sS(l-a;3-2^) 

X ja;^ + 22;B[-3a;^ + 2a;^x - 4a;x + 8a; - 3 - zx^^ ^ {x ^ x] 



„sc,yy _ ~vv _ { _z{l-x){x-%z-A) ^ a;(2x - 1) + 4^(x2 + 1) + 8^^ 

1 



Ci34 — ^123 — — h 



(1 — x)^a; (1 — a;)(l — a;g)a; 



2 



+ -— — \x^ + aj^a; — 6a; — 5a; + 3a;a; + 5 

(1 - x){x^ - Az) ^ 

+ z{x^ - 32x - 25x + 10a;x + 48) - 8^^(a; + x - 6)] 

32; 

+ 7 7T^^ rT^[82;^ - ^x^x - 20a; - 16x + 20a;x + 12 

(1 - a;)(a;2 - Azf 

+ Az{x^ - 10a; - 5x + 2a;x + 12) + 16^^] 

3x(l - a;) x{x^ + a; + 1) + 2z{^x^ - 2a; - 2) - '&z'^{x + 1) 



[x^ — AzY (1 — a;)a;- 



1 



^234 — Ci34 



3 r ^ /'^ _ ;?^2/'^2 I ~2 , I O;? _ _ 



_ zsS^ fa;p(a;-x)2(a;^ + x^ + 2a; + 2x-3)(l -a;) 



•^234 



I \ 1- Xg-AZ 



^ -[x^- 8a;^ + lla;^ + 90a;^ - 442a;^ + 908a;^ - 1008a;^ + 592a; - 72 



(1 Xg)Xg 

+ 2xx{2x^ - Ax^ - 80x^ + 477x^ - 1188a;^ + 1562a; - 532) 

+ x^x'^{imx^ - 771x^ + 2428a; - 1942) + Ax^x\-5x^ + 67x - 148) - 17a;%^] 

^-[x^ - 13x^ + 56x^ - llAx^ + 112a; - 22 

1 -a;p 

+ a;x(7a;^ - 58a;2 + 176a; - 114) + x'^x^{16x - 47)] 



28 



+ 7-^ h ^ a;) ^, 

1 -Xg J 

,sc,yy 1 — a: — 22; 2;(1 + x — Az) — x) 

+ 



1-x 2{l-xy 
1 + X + 2z(l -x) + z^{x - 2) + 2^3 



2(l-x + ^)(l-^)(l-x) 

^ -[{1 - x){l - X + X - 3xx) + z{-2x'^ 



2(1 - - z){l - Az){x^ - Az) 

+ llxx{-x + 2) - 18a; - 20x + 25) + z'^{x^ + 2a;x(2a; - 23) + 84a; + 68x - 60) 

+ 2z^{2xx{x + 3) + 16x + 21x - 124) + 4z\4xx - 24a; - 22x + 56) + 64/] 
3z 

+ 7i — r^[(l - x){-3xx + 2x + 5x-2) 

[l-x){x^-Azy 

xx{x — 12) + 28a; + 16x - 24) + Az'^{2x + x - 6)] 

B[(l - a;)(l + 4^(a: - 3)) + 2^^(2x^ - Ux + 11) + 8/(2a: - 1)] 

^ 2(1 - Az) ' 

7sc,yy _ isc,vv 

,sc,vv _ B{x-xy{l-Xg) _ 3x^[x'^ + xx-2{l-Xg)Y 
~ A{l-Xg-Az) 16(1 -x)(a;2 -4^)2 

tB 

[a;^ - lOx^ + 30a;^ - 38a;' - 6^2 



8(a;2 - Az) 

+ x^x^{Ax - 17) + xx{x^ + Ax^ - 29a;' + 54a; + 22x - 48) + 26a; + 12x - 6] 

+ — [2a;^ - 3a;^ - 5a;^ + 5a;'x' + 20a;' + a;x(2a;^ + 2a;' - 47a; + 40) - 64a; + 16 
16 

+ 8^S(-a;'' + 5a;^ + a;'x'(a; - 3) + 14a;' + a;x(a;^ - 4a;' - 3a; + 24) - 56a; + 22) 
+ 16z'S'(a;^ - 18x^ + x'x'(4a; - 3) + 62a;' + 6a;x(x' - 7a; + 9) - 88x + 24) 
+ 128a;g/5'(l - a;)'] + (a; ^ x), 

- ^"K^'V.^':^^^^ + . . J 2a;^ - 2a;^ - lla;^ 
8(1 - x){x^ - Azf 8(a;2 - Az) ^ 

+ a;'x'(6a; - 25) + 40a;3 - 52a;' - 8x' + 2a;x(a;^ + 2x^ - 20a;' + 38a; + 16x - 34) 
+ 36a; + 16x - 8] - ^[2a;^ + 2a;^x^ - 13a;^ + 22a;^ + 2a;'x'(3a;' - 5a; - 15) 



+ 8x'' - 40x^ + xx(6a;^ - 25a;^ - 8x^ + 120a; - 48) + 

+ 8zBx^g{x^x^ - a;^ + 12a;' + a;x(a;' - 11a; + 15) - 26x + 9) + IGz'^BXg] + {x ^ x), 



~yy ^ B{1 xg){x xf _B A _ 
'^ 4(1 -Xg- Az) 2x1 

+ 2a;x(2a;' - 12a; + 11) - 28a; + 6 + 2^a;'] + (a; ^ x), 
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_ Blx"^ + XX - 2{1 - Xg)]'^ B[2x^ - Sx"^ + xx{x - A) + 12x + Ax - 7] 
2s{2 - x){x^ - Az) ^ 2s{2 - x){l - X + z) 

- — [2a;V(4x - 11) + llx^ - 35x^ + a;x(-19x^ + 71x - 45) + 38a; - 7 

- Az{2x^x'^ - Ax^ + 16x^ + xx(3x^ - 17x + 17) - 23x + 6) 

- 82:^x^(1 — x)^] + (x ^ x), 

B{l-Xg-2z)xl\n{uj) Bx[x^ + XX - 2{1 - Xg)]'^ 
2/3(1 - ,Xg) ^ 8(a;2 - 4^) 

^2 

+ — [S(3x%^ + x^ + 2x3x3(2x^ - 33x + 78) - 20x^ + lOlx^ 

SXg 

+ x^x^(x^ - 28x^ + 215x2 - 668x + 535) - 226x^ + 272x2 

+ 2xx(-x^ + 22x^ - 128x^ + 322x2 - 431x + 152) - 176x + 24) 

- 4^Xg(3x^ - 7x2 ^ ^-^5^ _ 8) + _ 1) _ 16^25x^(2x2x2 - 4x^ + 16x2 

+ xx(3x2 - 17x + 17) - 23x + 6) - 32z^Sx2(l - x)2] + {x ^ x). (A.29) 



~AjA ~v^v^ 2s^^1 q ^'^^^g\ r 3 3 2 2 

7r(l - X)(l - Xg)(l - Xg - 42;) I 

2z 

+ xx(3x + 3x - 22) + 20x + 24x - 12 -\-x^ + 7x^ + 21x2 

1 — X 

+ xx(-x - 2x + 24) - 28x - 46x + 26 + 12^Xp]|, 

2[x + X - 3 + ^(x - 3) - ^2] + zx(4x + 1) + 12^2 



J/^^ _ AAA 



^123 ~ ^123 



1 — X (1 — x)(l — Xp)x^ 

^ (1 - 2x + Az){x^ - zx{%x + 1) + 8z2(6x - 1) - 16z=^(x + 5) + 64^^] 



2 



(1 - x)(l -Xg- Az){x - Az) 
[x^ - 15x^ + x^' + 2z(-2x^ + 52x=' + 3x2) 



(1 — x)(x — AzY^^ 

^2/ ■\ 2 .X 3/ 2 o 2Z^{\- X)\X^ \\2X - Azix \'0)\\ 

+ 8^2 _3™3 _ 28a;2 _ ^ 32z=^ x2 + 3x + 3 ^ — ^ ^ , 



%AA _ ~AA 
<-134 — <-123 

<-124 — <-124 



+ 2zss|(l - Xg)(x2 + XX - lOx + 5) - 2zS[x2x2 - 4x^ + 26x2 

+ xx(x2 - 22x + 31) - 44x + 11] - <oz^Bx\ + (x <^ x)|, 

sc,AA _ sc,yy , ~AA ~vv , o^J ^^(l + 12) , x^ + 2:x(4x + 1) + 12^2 

I 1 — xra^ 1 — X 1 — Xo x^ 
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+ xx{\ - x) - 2x + 1 + 1z{Zx^ - 9x + 2x + 5) + 1z^{-2x + x + 3) 

(1 -4z) 

+ 7:; ^^^^ + 2^^(2a^ - 5) + 2a; + 3x - 1 + ^(-5x^ 

{\ - x){x^ - ^zf 

+ a;x(a; - 10) + 20,t + 12,t - 16) + 4^^(2x + ,t - 5)] 

^ — 15a;^ + a; + 2z{x^ + 6a; + 6)] ^ a; + S — 3 + 2;(2x + x — 6) + 



(1 — x)x^ 1 — a; 

^234 — C134 _, 

41 = 44 + ^^^^1 - Y^[2x^ + x^x\-2x - 13) - 6a;5 - 12x^ 

+ x^x^{x^ - 22a;2 + 176,x - 250) + 98a;^ - 194a;2 + xx{x^ - 7x^ + 54a;^ 

- 258x^ + 542a; - 246) + 160a; - 24 - 4:z{3x^x^ - 4x^ + 42a;^ 

+ x^x'^{4x'^ - 66x + 144) - U8x^ + 262x'^ + xx{x'^ - 34x^ + IQOa;^ - 476a; + 274) 

- 232a; + 40) - 48^2x^(1 - a;)'] + ~ ""^^ [2x^ + x''x^{3x - 3) 

- 6a;^ + 12a;^ + a;x(a;^ - 3x^ - 6,t + 11) - 14a; + 3] + (a; ^ x)|, 

,sc,AA _ u^c,vv [ 3-^ (1 - .t)(1 - 2:) - 2:^ 2:(l + .x-42:) 22:^(1 -a;) 

Xl-x {l-x){l-z){l-x + z) {1-xf {1-xf 
B[{1 - x){-5 + 22z) + 2z'^{12x - 11) - 8^^] 
1-4^ 

1 , 2 - 

-[6a; — 2a;a; — 10a; + 5a; + 4 



(1 -a;)(l -^)(1 -4z)(x2 -4^)' 

- 2^(15a;^ - 5a;x - 23a; + 9x + 11) + 2z'^{12x^ + Axx - 2x - 15x + 4) 

6 

- 4:z^(4:xx + 8x- lOx - 5) - 16/1 - — rfa;x(-a; + 2) - x 

{l-x){x^-4zy 

- z{-6x'^ + xx{3x - 14) + 16a; + 12x - 10) 
+ 2/(a;^ + 2a;x - 12a; - 6x + 16) + 8z^]|, 



isc.AA 1 sc.AA 

^24 = ^13 



sc,AA ^ ,sc,yy z \ ^B{x - xf 3a;^[x^ + xx - 2(1 - Xg)f 

'14 0i4 ^ 4^\i-Xg-Az (l-x){x^-AzY 
2B 
x"^ — Az 



+ — ^[5a;^ - 21a;^ + 3x'^x^{x - 3) + 40a;3 - 38a;2 - 2x^ 



+ 2xx{Ax^ - 15x^ + 21x + 4x - 13) + 18a; + 4x - 2] - B^[B{x^x^{13x - 44) 

+ 7a;^ - 32x^ + x^x'^{7x^ - 58x^ + 128a; + 20) - 9a;^ + 150a;^ 

+ a;x(-14a;^ + 77a;^ - 44a;2 - 267a; + 214) - 180a; + 32) - 8z{x^ - 21x^ 
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xx{x + 10) + 56a; - 25) + 8z'^Bxg{x + x + 10){1 - xf] + {x ^ x)^ 

l,sc,AA _ usc,VV I Z f —3x[x'^ + XX - 2{1 - Xg)]'^ B ^^^4 ^^^3 



+ 4a;^x^ + 60x^ + 4,t^ + xx{12x^ - A2x - 9x + 52) - 52a; - 2Ax + 20] 
+ 25^ [3a; V -3x^ + 23x2 

+ xx(3a;2 - 29x + 33) - 40x + 10 + 4^2;^] + (x ^ x)|, 
rsc.AA _ rsc,vv 2zB{x — xY 

^sc,AA ^ ^sc,yy + — / - ^[^^ +^^-2(1 -Xg)]^ 2(l-x)(x^-xx-x + 2x-l) 
2s I (x2-42)(2-x) (1 -a; + ^)(2 -x) 

- B[B{x^x^ + x^ + x2x^(x2 - 16x + 28) - 12a;^ + 27a;^ 
+ a;x(-2a;^ + 26a;2 - 68a; + 33) - 22a; + 3) 

+ 4:z{7x'^ + Axx - 19x + 8) + 24^^Sa;g(l - a;)^] + (a; ^ x)|, 

.scAA ^ .scyy zf 2B{l-Xg- 2z)xl Injcv) B[x' + xx - 2(1 - x,)]^ 
2 1 P{l-Xg) (x^-Az) 

- B'^[B{-10x^x^ - lOx^ + x'^x^{-10x^ + 104a; - 149) + 64a;=^ - 120a;2 
+ 2a;x(10a;^ - 74a;2 + 162a; - 71) + 88a; - 11) + 2z(3a;^ - 40a;^ 

+ a;x(9a; - 34) + 106a; - 44) - 24:z'^Bxg{l - xf] + {x ^ x)|. (A.30) 
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Figure 4: The three jet fraction /| at NLO as defined in the text at = = mz as 
a function of Umin = Smin/ {si/cut) for the JADE algorithm at a value of the jet resolution 
parameter ycut = 0.03 with m^^(/i = mz) = 3 GeV and ^^(/x = mz) = 0.118. 
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Figure 5: Same as Fig. 0, but for the Durham algorithm. 
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Figure 6: The three jet fraction /g as defined in the text at ^/s = mz as a function of 
Ucut for the JADE algorithm. The dashed fine is the LO result. The NLO results are for 
^ = mz (solid line), ii = mz/2 (dotted line), and /i = 2m z (dash-dotted line). 
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Figure 8: Dependence of the three jet fraction /| for the JADE algorithm on the renor- 
mahzation scale at ^/s = mz and i/cut = 0.2 x 10^^^^° ^0.1 for on-shell b quark masses 
^poie _ 2 (j^ii curves) and ml"^^ = 5 GeV (dashed curves). 



42 



o 

1 



0.9 



0.8 



0.7 

50 100 150 



Figure 10: The double ratio C as a function of the cm. energy at ycut = 0.08 for the JADE 
algorithm. The full curve and the points show the NLO result, the dashed curve shows the 
LO result. In both cases a running mass m^^(// = y/s) evolved from m^^(/i = niz) — 3 
GeV is used. The dash-dotted curve shows the LO result for a fixed mass = 4.7 GeV. 
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Figure 11: The double ratio V as a function of i/cut for ^/s — /i — mz- Pull circles: NLO 
results for m^^(/i = mz) = 3 GeV. For comparison, the squares (triangles) are the LO 
results for rrib = 3 GeV {mi, = 5 GeV). The horizontal bars show the size of the bins in ycut- 
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Figure 12: Same as Fig. 11, but for the Durham algorithm. 
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Figure 13: Coefficient Ci of the thrust distribution (|6.71) for b quarks at y/s = mz (upper 
curve) and y/s = 30 GeV (lower curve). 
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Figure 14: Coefficient C2 of the thrust distribution (|6.71) for b quarks at y/s = mz (upper 
curve) and y/s = 30 GeV (lower curve). 
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